Landslides, threshold slopes and the survival of relict terrain in the wake of the Mendocino Triple Junction

Establishing landscape response to uplift is critical for interpreting sediment fluxes, hazard potential, and topographic evolution. We assess how landslides shape terrain in response to a wave of uplift traversing the northern California Coast Ranges (United States) in the wake of the Mendocino Triple Junction. We extracted knickpoints, landslide erosion rates, and topographic metrics across the region modified by Mendocino Triple Junction migration. Landslide erosion rates mapped from aerial imagery are consistent with modeled uplift and exhumation, while hillslope gradient is invariant across the region, suggesting that landslides accommodate uplift, as predicted by the threshold slope model. Landslides are concentrated along steepened channel reaches downstream of knickpoints generated by base-level fall at channel outlets, and limit slope angles and relief. We find evidence that landslide-derived coarse sediment delivery may suppress catchment-wide channel incision and landscape denudation over the time required for the uplift wave to traverse the region. We conclude that a landslide cover effect may provide a mechanism for the survival of relict terrain and orogenic relief in the northern Californian Coast Ranges and elsewhere over millennial time scales

Beyond the angle of repose: A review and synthesis of landslide processes in response to rapid uplift, Eel River, Northern California

In mountainous settings, increases in rock uplift are often followed by a commensurate uptick in denudation as rivers incise and steepen hillslopes, making them increasingly prone to landsliding as slope angles approach a limiting value. For decades, the threshold slope model has been invoked to account for landslide-driven increases in sediment flux that limit topographic relief, but the manner by which slope failures organize themselves spatially and temporally in order for erosion to keep pace with rock uplift has not been well documented. Here, we review past work and present new findings from remote sensing, cosmogenic radionuclides, suspended sediment records, and airborne lidar data, to decipher patterns of landslide activity and geomorphic processes related to rapid uplift along the northward-migrating Mendocino Triple Junction in Northern California. From historical air photos and airborne lidar, we estimated the velocity and sediment flux associated with active, slow-moving landslides (or earthflows) in the mélange- and argillite-dominated Eel River watershed using the downslope displacement of surface markers such as trees and shrubs. Although active landslides that directly convey sediment into the channel network account for only 7% of the landscape surface, their sediment flux amounts to more than 50% of the suspended load recorded at downstream sediment gaging stations. These active slides tend to exhibit seasonal variations in velocity as satellite-based interferometry has demonstrated that rapid acceleration commences within 1 to 2 months of the onset of autumn rainfall events before slower deceleration ensues in the spring and summer months. Curiously, this seasonal velocity pattern does not appear to vary with landslide size, suggesting that complex hydrologic–mechanical feedbacks (rather than 1-D pore pressure diffusion) may govern slide dynamics. A new analysis of 14 yrs of discharge and sediment concentration data for the Eel River indicates that the characteristic mid-winter timing of earthflow acceleration corresponds with increased suspended concentration values, suggesting that the seasonal onset of landslide motion each year may be reflected in the export of sediments to the continental margin. The vast majority of active slides exhibit gullied surfaces and the gully networks, which are also seasonally active, may facilitate sediment export although the proportion of material produced by this pathway is poorly known. Along Kekawaka Creek, a prominent tributary to the Eel River, new analyses of catchment-averaged erosion rates derived from cosmogenic radionuclides reveal rapid erosion (0.76 mm/yr) below a prominent knickpoint and slower erosion (0.29 mm/yr) upstream. Such knickpoints are frequently observed in Eel tributaries and are usually comprised of massive (> 10 m) interlocking resistant boulders that likely persist in the landscape for long periods of time (> 105 yr). Upstream of these knickpoints, active landslides tend to be less frequent and average slope angles are slightly gentler than in downstream areas, which indicates that landslide density and average slope angle appear to increase with erosion rate. Lastly, we synthesize evidence for the role of large, catastrophic landslides in regulating sediment flux and landscape form. The emergence of resistant blocks within the mélange bedrock has promoted large catastrophic slides that have dammed the Eel River and perhaps generated outburst events in the past. The frequency and impact of these landslide dams likely depend on the spatial and size distributions of resistant blocks relative to the width and drainage area of adjacent valley networks. Overall, our findings demonstrate that landslides within the Eel River catchment do not occur randomly, but instead exhibit spatial and temporal patterns related to baselevel lowering, climate forcing, and lithologic variations. Combined with recent landscape evolution models that incorporate landslides, these results provide predictive capability for estimating erosion rates and managing hazards in mountainous regions

Oregon 2100: projected climatic and ecological changes

Greenhouse climatic warming is underway and exacerbated by human activities. Future outcomes of these processes can be projected using computer models checked against climatic changes during comparable past atmospheric compositions. This study gives concise quantitative predictions for future climate, landscapes, soils, vegetation, and marine and terrestrial animals of Oregon. Fossil fuel burning and other human activities by the year 2100 are projected to yield atmospheric CO2 levels of about 600-850 ppm (SRES A1B and B1), well above current levels of 400 ppm and preindustrial levels of 280 ppm. Such a greenhouse climate was last recorded in Oregon during the middle Miocene, some 16 million years ago. Oregon’s future may be guided by fossil records of the middle Miocene, as well as ongoing studies on the environmental tolerances of Oregon plants and animals, and experiments on the biological effects of global warming. As carbon dioxide levels increase, Oregon’s climate will move toward warm temperate, humid in the west and semiarid to subhumid to the east, with increased summer and winter drought in the west. Western Oregon lowlands will become less suitable for temperate fruits and nuts and Pinot Noir grapes, but its hills will remain a productive softwood forest resource. Improved pasture and winter wheat crops will become more widespread in eastern Oregon. Tsunamis and stronger storms will exacerbate marine erosion along the Oregon Coast, with significant damage to coastal properties and cultural resources

Frost for the trees: Did climate increase erosion in unglaciated landscapes during the late Pleistocene?

Understanding climatic influences on the rates and mechanisms of landscape erosion is an unresolved problem in Earth science that is important for quantifying soil formation rates, sediment and solute fluxes to oceans, and atmospheric CO2 regulation by silicate weathering. Glaciated landscapes record the erosional legacy of glacial intervals through moraine deposits and U-shaped valleys, whereas more widespread unglaciated hillslopes and rivers lack obvious climate signatures, hampering mechanistic theory for how climate sets fluxes and form. Today, periglacial processes in high-elevation settings promote vigorous bedrock-to-regolith conversion and regolith transport, but the extent to which frost processes shaped vast swaths of low- to moderate-elevation terrain during past climate regimes is not well established. By combining a mechanistic frost weathering model with a regional Last Glacial Maximum (LGM) climate reconstruction derived from a paleo-Earth System Model, paleovegetation data, and a paleoerosion archive, we propose that frost-driven sediment production was pervasive during the LGM in our unglaciated Pacific Northwest study site, coincident with a 2.5 times increase in erosion relative to modern rates. Our findings provide a novel framework to quantify how climate modulates sediment production over glacial-interglacial cycles in mid-latitude unglaciated terrain

A 1-D Mechanistic Model for the Evolution of Earthflow-Prone Hillslopes

In mountainous terrain, deep‐seated landslides transport large volumes of material on hillslopes, exerting a dominant control on erosion rates and landscape form. Here, we develop a mathematical landscape evolution model to explore interactions between deep‐seated earthflows, soil creep, and gully processes at the drainage basin scale over geomorphically relevant (\u3e103 year) timescales. In the model, sediment flux or incision laws for these three geomorphic processes combine to determine the morphology of actively uplifting and eroding steady state topographic profiles. We apply the model to three sites, one in the Gabilan Mesa, California, with no earthflow activity, and two along the Eel River, California, with different lithologies and varying levels of historic earthflow activity. Representative topographic profiles from these sites are consistent with model predictions in which the magnitude of a dimensionless earthflow number, based on a non‐Newtonian flow rheology, reflects the magnitude of recent earthflow activity on the different hillslopes. The model accurately predicts the behavior of earthflow collection and transport zones observed in the field and estimates long‐term average sediment fluxes that are due to earthflows, in agreement with historical rates at our field sites. Finally, our model predicts that steady state hillslope relief in earthflow‐prone terrain increases nonlinearly with the tectonic uplift rate, suggesting that the mean hillslope angle may record uplift rate in earthflow‐prone landscapes even at high uplift rates, where threshold slope processes normally limit further topographic development

Climatic controls on frost cracking and implications for the evolution of bedrock landscapes

of bedrock landscape

Sediment yield, spatial characteristics, and the long-term evolution of active earthflows determined from airborne LiDAR and historical aerial photographs, Eel River, California

In mountainous landscapes with weak, fine-grained rocks, earthflows can dominate erosion and landscape evolution by supplying sediment to channels and controlling hillslope morphology. To estimate the contribution of earthflows to regional sediment budgets and identify patterns of landslide activity, earthflow movement needs to be quantified over significant spatial and temporal scales. Presently, there is a paucity of data that can be used to predict earthflow behavior beyond the seasonal scale or over spatially extensive study areas. Across 226 km^2 of rapidly eroding Franciscan Complex rocks of the Eel River catchment, northern California, we used a combination of LiDAR (light detection and ranging) and orthorectified historical aerial photographs to objectively map earthflow movement between 1944 and 2006. By tracking the displacement of trees growing on earthflow surfaces, we find that 7.3% of the study area experienced movement over this 62 yr interval, preferentially in sheared argillaceous lithology. This movement is distributed across 122 earthflow features that have intricate, elongate planform shapes, a preferred south-southwesterly aspect, and a mean longitudinal slope of 31%. The distribution of mapped earthflow areas is well-approximated by a lognormal distribution with a median size of 36,500 m^2. Approximately 6% of the study area is composed of earthflows that connect to major channels; these flows generated an average sediment yield of 19,000 t km^(−2) yr^(−1) (rock erosion rate of ∼7.6 mm/yr) over the 62 yr study period, equating to a regional yield of 1100 t km^(−2) yr^(−1) (∼0.45 mm/yr) if distributed across the study area. As such, a small fraction of the landscape can account for half of the regional denudation rate estimated from suspended sediment records (2200 t km^(−2) yr^(−1) or ∼0.9 mm/yr). We propose a conceptual model for long-term earthflow evolution wherein earthflows experience intermittent activity and long periods of dormancy when limited by the availability of readily mobilized sediment on upper slopes. Ultimately, high-order river channels and ephemeral gully networks may serve to destabilize hillslopes, controlling the evolution of earthflow-prone terrain

Raw elevation grids of experimental eroding landscape

The zeroth -order response of eroding topography to changes in base-level is driven by patchy and intermittent mass flux, both on hillslopes and in channels. However, most models of landscape evolution use continuously differentiable equations that average over these patchy first-order transport processes. Because of the limited time- and space-resolution of field observations, the relationship between zeroth-order landscape evolution and first-order fluctuations about the zeroth-order state due to patchy, intermittent transport remains unclear. Here, we use five physical experiments of an eroding experimental landscape to examine how the signature of first-order transport, as described by autocorrelation functions of local elevation time series, varies as a function of the vigor of hillslope transport relative to channel incision. Our results show that experiments with higher hillslope transport efficacy have higher autocorrelation coefficients, suggesting that differences in zeroth-order transport coefficients may be driven by differences in patchy, first-order transport processes. These higher autocorrelation coefficients also imply that in landscapes where hillslope transport dominates, landscape dynamism is reduced and landforms are more persistent over time